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Abstract 

We show that the lancu-Mueller factorization has a simple interpretation in the 
Reggeon - hke technique based on the BFKL Pomeron. The formula for calculating 
the high energy asymptotic behaviour for the colour dipole-dipole amplitude is pro- 
posed which suggests a procedure to calculate this amplitude through the solution to 
the Balitsky-Kovchegov non-linear equation. We confirm the lancu - Mueller result that 
a specific set of enhanced diagrams is responsible for the high energy behaviour for fixed 
QCD coupling. However, it is argued that in the case of running QCD coupling, this 
asymptotic behaviour originates from the Balitsky-Kovchegov non-linear equation. A 
new solution to the non-linear equation is found which leads to a different asymptotic 
behaviour of the scattering amplitude even for fixed as- 
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1 Introduction 

lancu and Mueller in recent papers [Ij have suggested a new approach^ to determine the 
high energy asymptotic behaviour of the scattering amplitude which goes beyond the non- 
linear equation (BK equation 011]) of the dense partonic system. In this system the gluon 
occupation numbers are large, the gluonic fields are strong and such a system enters a new 
phase of QCD: the colour glass condensate 0011111! (see also Ref. for a recent review of this 
approach). The key new element in lancu-Mueller approach is an attempt to take into account 
fluctuations in the partonic wave function of the fast moving particle which were neglected 
in the non-linear equation. These fluctuations could be taken into account by more general 
approach to CGC related to so called JIMWLK equation, which is functional equation and 
it is not very practical at the moment. We believe that the lancu-Mueller approach gives as an 
opportunity to find more transparent and analytic way to take into account the fluctuations 
in the partonic wave function. 

The main goal of this paper is to show that the lancu-Mueller approach arises naturally 
from the Reggeon - like diagram technique [Hj based on the BFKL Pomerons pH]- This 
technique is the well known way of incorporating the fluctuation in the partronic wave function 
for high energy scattering which leads to the non-linear BK evolution equation PT|. We will 
review the region of applicability of the non-linear equation for the scattering processes. It 
was discussed in Ref. j3] (see also Refs. fT^ IT^ [Hj) but it has not been utilized since that 
time. We show that the lancu-Mueller result for the scattering amplitude is valid for frozen 
QCD coupling in a large but limited range of energies, while for running QCD coupling we 
expect rather the answer obtained in the colour glass condensate. 

We concentrate our efforts on understanding dipole-dipole scattering. The main idea of 
the lancu-Mueller factorization as well as our approach to this factorization can be illustrated, 
considering the first so called enhanced diagram (see Fig. ^a) which describes the fluctuation 
in the partonic wave function of the incoming fast dipole. 

We calculate this diagram in the simple toy - model for the BFKL ladder, namely, assuming 
for the BFKL Pomeron with the intercept A, a simple expression e^^, while neglecting any 
dependence on the size of the interacting dipoles jTHj. In addition, we impose the relation 
Gi = gi = A and G2 = 92 = ^ where Nc is the number of colours. These relations follow 
directly from the QCD estimates [13 HE! , but the most important Nc suppression is the old 
result of topological expansion [T7j . 

The expression for Fig. ^a has the form: 



^The key ingredients of this approach have been suggested by Mueher and Salam 2 but at that time the 
approach based mostly on numerical simulations while now we are able to develop analytical methods. 
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Fig. [Ba Fig. [T}b 

Figure 1: The first enhanced diagram (Fig. [T]-a) for dipole-dipole scattering and the lancu- 
Mueller factorization for this diagram (Fig. ^b) 



In Eq. (jl.ip the sign minus reflects the fact that such diagrams describe shadowing which 
tames the power-hke energy increase e^^, caused by the BFKL Pomeron exchange. For large 
Nc there exists a region of energy given by the inequality 

N'. » ^e^^ > 1 (1.2) 

in which only the first term in Eq. (jl.lj) is large (of the order of unity), while all other terms 
are small. In this region of energy the diagram has a simple form 

A{Fzg.m-a) = (-1)^6^^^ (1-3) 

c 

The main idea of lancu and Mueller [T] is to calculate the amplitude of Eq. ()1.3p using a 
different approach, which is presented graphically in Fig. ^b. This approach is based on the 
generating function [13 HEl HHj 

Z{y,u) = Pn{Y-y)u'' (1.4) 

n 

where Pn{y) is the probability of having n-dipoles inside the fast dipole. We will show in 
section 3 that in the simple model, in which we have calculated the diagram of Fig. ^a, the 
term in the generating function which is responsible for the amplitude for production of two 
dipoles (see Fig. [TJ-b) has the following form: 

A{Y-y,^) = 1- Z{Y -y,u = 1 + -/) = (-1) 7' e^^^^'^^ (^e^iy-v) _ 1) (1.5) 
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Therefore, Eq. (jl.3p can be written in a different form, namely 
AiFig.m-a) = (-1) A(r - |/, 7 = 1) 7 = ^) 



A' {-m{Y - y) P,{y) (1-6) 
(-1) A2 d' (1 -Z{Y-y,l + 7)) , (1 - Z{y, 1 + j/N^)) 



=0 717^. 17= 



4 (d7)2 (d7)^ 

Eq. ()1.6|) leads to the answer, which is not the same as the correct expression for the enhanced 
diagram given in Eq. (jLip . but it correctly reproduces the leading term (see Eq. ()1.3p ). The 
corrections to this term depend on the choice of the value for the rapidity y (see Fig. E-b) 
which has only auxiliary meaning. The fact that the generating function Z depends on y 
is very transparent, since the wave function of the dipole depends on the reference frame. 
However, the amphtude as a physical observable should not depend on y. It does not depend 
on y provided 

^4 e - « 1 and e ~ « 1 (1.7) 

The choice y = Y/2 leads to the estimates that both terms in Eq. ()1.7p are of the order of 
-C 1, and is most accurate. This simple seems to justify the choice of the reference 
frame suggested in Ref. [1. 

Formula of Eq. p.6|) is the simplest example of the lancu- Mueller factorization. In section 
3 we will consider in detail a toy-model which we have used here to illustrate the main idea 
of the lancu-Mueller factorization. We will show that Eq. ()1.6|) has a natural generalization 
which coincides with the direct sum of all enhanced diagrams in the kinematic region given 
by Eq. (jl.2j) . In section 4 the model result will be generalized for the QCD case to include a 
dependence on the sizes of the interacting dipoles. 



2 Non-linear equation 

The non-linear equation the final form of which at fixed impact parameter was suggested 
by Bahtsky [S] and Kovchegov sums all so called 'fan' diagrams (see Fig. |21-a). We would 
like to repeat here the arguments that led to this equation EI for the case of the dipole-dipole 
scattering. The derivation given in Refs. [Sill] used the fact that the target is a heavy nucleus. 
As it will be clarified later, the nuclear target indeed enlarges the region of applicability for 
the non-linear equation, but it is not essential for the saturation in the colour glass condensate 
domain. 



2.1 The first enhanced diagram: 

For a deeper understanding of the non-linear equation we start with the calculation of the 
same enhanced diagram of Fig. ^a. The gluon 'ladder' shown in this figure is the solution to 
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Figure 2: The first 'fan' diagram (Fig. |21-a) for dipole-dipole scattering that is summed by 
the non-hnear evolution equation, and the first diagram that is not included in summation by 
the non-linear equation (Fig. |21-b) 

the BFKL equation [TU], which can be written for the dipole (ri)-dipole(r2) amplitude in the 
factorized form HH Cg 1201 EH (see Fig. H ). 



at 



(2.8) 



27ri 



D{u) e-(^)^V^(ri,i,g;i/)V^(r2,t,g;-z/) 



with 



and with 



asNc 



TT 



x{i) 



D{v) = 

asN, 



,2 I 1^2 



TT 



(2^(1) - V^(7) - z^(l-7)) 



(2.9) 
(2.10) 



where ip{f) = dlnT{f)/df , r(/) is the Euler gamma function, 7 = ^ — iu and where 

27r2 



r? 



{R^ + hri,ty {Ri - ^n,ty 



(2.11) 



and the following is our notation: y = ln(xo/x); rj.f is the size of the colour dipole "z" and 
Ri is the position of the center of mass of this dipole. q is momentum transferred along the 
BFKL Pomeron (see Fig. El). 

The coefficient 6(z/) in Eq. ()2.11|1 is defined in Refs. jTHlEII, but we do not need its explicit 
form in what we discuss below. The only important property of Eq. (j2.1Hl that we will use is 



5 



N(ri,r2;y,q) = 




= q 




V{ri ,q ; V ) 
co(v)y 



V(r2,q;^) 



Figure 3: The BFKL Pomeron. 



[mini 



V {rt,q;i^) ,V {rt,q;-iy,) = 6{iy - u') 



(2.12) 



The simplest enhanced diagram of Fig. 0] has been given in Ref . j2I] and it is equal to 



8 V 27r2 



j (f'k d'^ d'^i J dyi J dy2 (2-13) 



V^(ri,i, q- u)D{u) e^^^^^ (^-^^^ Gj,p{y- z/i, v^, q, k) D{u^) Diu^) e(-(-^) + -('^^))(J'^-^^)G3p(z/; Ui, 1^2] q, k) 

where 7, = 1/2 + iz/j and the triple Pomeron vertex G^p is calculated in Ref. [21] all details 
about Eq. ()2.13p as well as its explicit derivation, is given. 

In spite of the additional complicated integrations over momenta transferred and anoma- 
lous dimensions 7^, appearing in Eq. (j2.13p . the equation has the same integrations over 
rapidities as in the simple toy model. Calculating the integrals over yi and y2 and keeping 
only the maximal power of energy, we reduce Eq. ()2.13p to a simpler expression 



^,4 



a 



32 m 



yi 



dkd-id-iid^2 / dyi / dy2 V{ri^t,q;iy) (2.14) 



Diiy) G3p{u; z/i, 1/3; q, k) D{vi) D{v2) 



>(.0+-(-2))yG'3p(zy;i.i,i.2;g,fc) 



D{u) V{n^uq-u) 

Comparing Eq. ()2.14|) with Eq. ()2.8|) one can see that the ratio of these diagrams is of the 
order of 

N{Fig.^ 



At large values of Y this ratio becomes of the order of unity, and to obtain a correct scattering 
amplitude the enhanced diagrams should be taken into account. 
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2.2 'Fan' diagrams: 



At first sight this discussion leaves no room for a special status of the 'fan' diagrams (see 
Fig. 121) that are summed by the non-linear equation jHlll]- Indeed, as was analyzed long ago 
[Sj, we cannot expect any suppression of the enhanced diagrams for the fixed QCD coupling, 
which has been used in the above calculations. However, the situation changes crucially if we 
consider the scattering amplitude for two dipoles with different sizes (say ri ^ = r -C r2^t = R 
) and we take into account a running QCD coupling. For rather rough estimates let us assume 
that we can include running ashy replacing Eq. ()2.10|) by 



(2.16) 



G3p(V;Vi,V2) 

JX)(A^)(yi-y2) 



G3p(Vi,V2;V) 




V(ri,q ;V) 

CO(v)(Y-yi) 



.co(Vi)(yi-y2) 



C0(v)y2 



V(r2;q,^) 



Figure 4: The first enhanced di- 
agram in the Reggeon calculus 
with the BFKL Pomerons (zigzag 
lines) . 



In this case, the typical value for Y — yi in Eq. ()2.13|) is of the order of 



Y - yi 



1 1 u-1 1 

^ 1 while ?/2 ~ 



as{ri, 



<^ Y - yi 



(2.17) 



If ^2,t ~ R where i? is a typical hadron size, as{R) ~ 1 and y2 ~ 1. For such small 
y we cannot trust the Reggeon like diagrams, but we can replace the enhanced diagram of 
Fig. ^or of Fig. E]by the 'fan' diagram (see Fig. |21). The meaning of such a replacement is 
that in a 'fan' diagram the interaction of BFKL ladders with the target, is a subject of the 
non-perturbative QCD approach. Assuming that the BFKL ladders interacts with the target 
independently (without correlations), we can write the non-linear equation. 

Therefore, the non-linear equation is a direct consequence of two physical ideas: (i) the 
high energy amplitude can be replace by exchange of the BFKL 'ladder'; and (ii) in the first 
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approximation the correlations can be neglected for the low energy interactions of the BFKL 
'ladder' with the target. The last claim has a theoretical justification for a nucleus target jl], 
here we do not consider the interaction with nucleus. 



2.3 Generating functional: 

In Ref. jl5^ A.H. Mueller suggested separating the high energy part of 'fan' diagrams from 
the low energy part, by introducing the generating functional. This functional has the form: 

/n 
Pn{Y ~ y,r;ri,bi;r2,b2; . . . n^bi] . . . r„,6„) J| u{fi,bi) (f n (fbi 

n=i i=l 

(2.18) 

where u{ri) = Ui is an arbitrary function of rj and bi. Pn denotes the probability density 
of finding n dipoles with rapidity y, with transverse size ri,r2, . . . rj . . . r„, and with impact 
parameters &i, &2, • • • h ... 6„ with respect to the mother-dipole, in the wave function of the 
fast moving dipole of the size r and rapidity Y > y. The functional of Eq. ()2.18p satisfies two 
conditions: 

1. At y = Y, Pi = 5^'^\f — fi) with Pn>i = 0. This condition means that at the 
beginning of the evolution we have one fast moving dipole or, in other words, we sum the 
'fan' diagrams which start with the exchange of one BFKL Pomeron. For the functional 
we have 

Z(r-y = 0,r;H) = «(r) ; (2.19) 

2. At Mi = 1 

ZiY-y,r;[u,= l]) = 1; (2.20) 

Eq. fl2.2()|l expresses the physical meaning of the functional:the sum over all probabilities 
is equal to unity. 

This functional sums all 'fan' diagrams (see Fig. El ), and the simple linear functional 
equation can be written for it 



dZiY-y,r,;[u,]) _ . n . x <^ 



das y 



J d\ u{n) uj{n) — Z{Y-y,r; [u^]) (2.21) 



where 



(fri (f r' u{ri) u{fi - -^—^ Z(Y -y,r; [u{r'),Ui]) 

u [Ti — ry ou[r) 

-C) = ^/<i^/3.7^. (2.22) 



the notation 5/5ui is used for the functional derivative. 

The physical meaning of each term is clear: the first one describes a probability for the 
BFKL Pomeron to propagate from rapidity y to rapidity y + dy without any decay, while the 
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Figure 5: The 'fan' 
diagrams that are 
summed by the gen- 
erating functional. 
Zigzag hues denote 
the BFKL Pomerons. 



second term accounts for the possibihty for decay of one dipole to two dipoles or, in other 
words, it takes into account the triple BFKL Pomeron vertices. 

A general property of Eq. ()2.2H) . is that a solution to this equation can be written as a 
function of a single variable Rewritmg dZ/dy a.s dZ/dy = J d?r {5Z/5u{r)) {du{r)/dy) 
and using the initial condition of Eq. ()2.19|) one obtains a non-linear equation ^3] : 

__ = -u[r) ZiJ -y,r\\u^\) (2.23) 

+ / ^^^ %.2(/lp)2 Z{y-Vy-An^s) Z(F-y,r-r-');M . 

Eq. ()2.21|) and Eq. ()2.23|) solve the problem of finding the probability to produce a number 
of dipoles of different sizes at rapidity from the single fast moving dipole at rapidity Y . 
This probability is independent of the target. Assuming that all produced dipoles interact 
with the target independently (without correlations), we can calculate the resulting scattering 
amplitude for the 'fan' diagram (see Fig. |2I). As was proved in Ref. this amplitude is equal 
to 

N{Y,r-Yl{rM\) = 1 " ^ (F, r, 6^; [7(r,,, 6,) + 1]) . (2.24) 

where —'y{ri,bi) is the amplitude of the interaction of the dipole of size rj, at impact parameter 
b and at y = with the target. 



2.4 The lancu-Mueller factorization and enhanced diagrams: 

The lancu-Mueller approach to calculations of the enhanced diagrams contribution to the 
scattering amplitude is shown in Fig. IHl We need therefore 
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(i) to calculate the amplitude of production of n-dipoles of the different sizes, that are 
produced by both colliding dipoles with sizes r and R in Fig. IHl 

(ii) to multiply the product of these amplitudes by 7(rj,r^) for each pair; 

(iii) to integrate over all r^, bi and r^, b[. 

Summing over all possible numbers of the interacting dipoles give us the amplitude. This 
procedure leads to the following formula (see Fig. IHl): 

oo „ n n 

Nir,R,Y;b) = J U d% K U d\ r[ ^{rj - h;r[,b[) (2.25) 

n=l i=l i=l 

Nn{r,Y -y;ri,b-bi;r2,b-b2; . . . rib -bf, . . . b -b^) Nn{R,y;r[,b\;r2,b2; . . . r[,b[ . . . r 
where Nn is the general term of expansion of the amplitude given by Eq. ()2.24|1 

n 

~ H n d^^id^h 'y{ri,bi)Nnir,Y;ri,bi;r2,b2 . . . ri,bi; . . . rn,bn) ■ (2.26) 

n=l i=l 

In Eq. ()2.26p we view the scattering amplitude N{Y,r] [7(rj,&j)]) as a generating functional 
with respect to arbitrary functions 7(rj). To obtain the scattering amplitude for the particular 
process we have to replace functions 7(rj) by 7(rj, R) where 7(rj, R) = —A{dipole — target) 
where A is the amplitude of the dipole - target interactions at low energies with sizes and 
R, respectively. 

Comparing Eq. (|2.24p and Eq. (j2.26j) one can see that Nn can be calculated as 

1 " (5 

Nn{r,Y,b,ri,bi;r2,b2; . . . ri,bi; . . . rn,bn) = —{Ti^C^ - Z {Y,r;[j{ri,bi) + 1]) ) \^^=o 

(2.27) 

As we have discussed, the artificial rapidity y cancels in the product of Eq. ()2.25p . 

In principle Eq. ()2.25p solves the problem of calculating the scattering amplitude, if we can 
find the expression for 7(rj,r^). Comparing Eq. ()2.8|) and Eq. ()2.13p and taking into account 
the completeness relation of Eq. ()2.12j) we obtain that 

m,K^r,,r'^ = (-l)'^^'n^.-^:)^^'Hr-;-r-:)^-^l (2.28) 

Returning to the impact parameter representation we have Eq. ()2.28j) in the form 

7(6.,6:;r„r:) = i-^) ^^'Hn - r^} ' ^ ^ (2.29) 

We recall that ds = Ncas/ir. 
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We can clarify the physical meaning of Eq. ()2.28|) and Eq. ()2.29|) by considering the first 
term in Eq. ()2.27|1 . which corresponds to the single BFKL Pomeron exchange. This term is 
equal which corresponds to the 

iV(F-y,r,g;[7(ri)]) = / rfViiV^^^^ (y - r, n; g) 7(ri) (2.30) 

Choosing 7(ri) = v and ^{r'x) = 7(q'; "^x-, f'l)/^^ one can see that Eq. ()2.2(jj) could be rewritten 
in the form 

^BFKL ^. ^) ^ _ j ^2^^ ^2^/ ^ - y, r, b- [7(ri)]) N {y, R, q; [-f{r[)]) (2.31) 
= - / d'n d\[N^^^^ {Y - r, n-q)u N^^^^ {y, R, r[- q) iMlillil 
= _ I rfVi d\[N^^^^ {¥ - r, n; q) N^^^^ (y, R, r[; q) ~^{q- n, r[) 

In the last equation we use the completeness of the BFKL vertex functions V{r, q, v) (see 
Eq. ()2.12|1 ) and the explicit form of Eq. ()2.28j) for 7(g; ri, r'^), to take the integrals over and 
r[. Eq. ()2.3H) gives the main formula that allows us to reduce the product of two functional to 
the exchange of the BFKL Pomeron and, therefore, to the Reggeon -like diagram technique. 

This formula together with Eq. ()2.25|) and Eq. ()2.27|) leads to the calculation of the scat- 
tering amplitude in the lancu-Mueller approach. However, before performing a calculation 
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using Eq. (|2.25p and Eq. (j2.27p we would like to clarify our approach using a simple model 
suggested in Ref. P^ . 



3 The lancu- Mueller factorization in a toy model 

We have used this model in the introduction to illustrate the main points of our approach to 
the lancu-Mueller factorization. In the toy model, the probability for the dipole to decay in 
two dipoles is a constant (ujq), which is equal to the probability for a dipole to survive without 
producing another dipole. In other words, in the framework of the Reggeon -type diagram 
technique, the intercept of Pomeron is equal to G^p - vertex, and both are equal to ujq. In 
this model we also neglect the fact that we have dipoles of different sizes, and therefore the 
generating functional of Eq. ()2.18|) reduces to a generating function 

Z{Y-y,u) = P^u'^ (3.32) 



n=l 



with two initial and boundary conditions: 



Aty = Y: Z{Y-y = {},u) = u ; (3.33) 
Atu=l: ZiY-y,u=l) = 1; (3.34) 



Eq. (jTTTIl reduces to (see Refs. [131111; 



dZ{y,u) dZ{y,u) 

- uJo{u{l-u)} (3.35) 



dy du 

which has the solution 



u 



Using Eq. (j2.24p one obtains 



^ pojo (Y-y) 

N(Y-y,-f) = (3.37) 



Eq. ()3.37|) together with a natural reduction of Eq. ()2.29|) to 

7 = ISM = (-1)4? (3-38) 

allows us to use Eq. ()2.25|1 to estimate the high energy asymptotic behaviour of the scattering 
amplitude, in this simple model. However, we suggest a more compact formula than Eq. ()2.25|1 . 
namely, 

NiY) = ^ <f —NiY-y,^sMi^)N{y,-); (3.39) 

2 TT « JCi V V 
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Figure 7: Complex v plane and 
contours for integration over v in 
Eq. and Eq. . 



where the integration contour is a circle with the unit radius around z/ = (see contour C\ 
in Fig. [7j). Expanding the functions Z with respect to v and \ jv one can prove in the case of 
our simple model that Eq. ()3.39p is equivalent to Eq. ()2.25p . 

Taking the integral of Eq. (|3.39|) explicitly we obtain the answer 



NiJ) = = — ^ (3.40) 

1 + ^ ( e'^oiY^y) - 1) ( e'^oy - 1) 

where ■jsm is given by Eq. (|3.38|) . On can see that Eq. (|3.4Up leads to 

limiV(y) ^ 1 - ^e--«^ + o( ^e"-oy/2\ (3 41) 

4 High energy asymptotic behaviour of the scattering 
amphtude 

In this section we return to the general expression for the lancu-Mueller factorization given by 
Eq. (j2.25j) and Eq. (j2.27j) . We would like to calculate the high energy asymptotic behaviour 
based on our experience with the simple model, that has been discussed in the previous section. 

We start with discussion of dipole-dipole amplitude in the simplified model for the BFKL 
kernel, namely, assuming that 

r i ioTr'Ql<l; 
^(^) = (4.42) 
[ ^ for r2g2 > 1- 
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instead of the general BFKL kernel given by Eq. (j2.1(Jj) . In Eq. (j4.42j) denotes the saturation 
momentum. It is shown in Ref . [22 that Eq. ()4.42|) sums double log contributions of the order 
of [asY \a{r'^ A^))^ in the kinematic region of perturbative QCD, namely, r"^ Q\ < 1, while in 
the saturation region (r^ > 1) it takes into account large terms of (05 ln(r^ Q^))" -type. 
In section 4.3 we will return to discussion of the general kernel of Eq. ()2.1()|1 . 

4.1 Fixed ols'. 

4.1.1 Generating functional for the Balitsky - Kovchegov scattering amplitude: 

We calculate the high energy asymptotic behaviour of the scattering amplitude using our 
experience with the simple toy model, and the fact that amplitude N {Y, r, b; [7(^4)]) has been 
found in Ref. [22j. We will show that a direct generalization of Eq. ()3.39|1 leads to the high 
energy amplitude that has been proposed in the lancu-Mueller paper P, namely, 

N{Y,ri,r2;b) =^ 1 - e'^^^^-^"^' (4.43) 

where coefficient c determines the asymptotic behaviour of the solution to the Balitsky- 
Kovchegov non-linear equation. This behaviour was found in Ref. [22] and it has the form 

N'^^{Y,n;b) =^ 1 - e-^(^-^o)' (4.44) 

with c = 2 ag. 

The first problem that we need to solve, is to build the generating functional for the 
amplitude (see Eq. ()2.26|) ). based on the solution given in Ref. [22] • 

n < l/Q4Y-y,b) 

As was shown in this paper, the solution for the short distances {vtQs < 1) is the single 
BFKL Pomeron, and, therefore, the generating function for the scattering amplitude at these 
distances is 

Nsd{Y -y,ri,b;[^{n)]) = j d\ ^{n) N''^''\{Y - y,n,n,b) (4.45) 

n > l/QsiY-y,b) 

For long distances, the solution for the interaction of the dipole (ri) with the dipole (r^) 
can be written in the form p2j: 

N{Y -y,ruri;b) = 1 - e''^^'^ (4.46) 
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where 

z = \n{rlQl) = Aas{Y-y) - Hr^/rf) + 2/5(6) (4.47) 

for the kernel of Eq. ()4.42|1 . j3{b) = 2 In S{b) in this approach, where S{b) is given by the 
impact parameter dependence of the Born approximation amphtude for this interaction (see 
Ref. [22] for details). 

Function is defined [22] by 

.^V-2f . (4.48) 

where the value of 0o should be found by matching the solution with Eq. (j4.46p at ri = l/Qs- 
Assuming that 0o is small we obtain from Eq. (j4.48p that for z <\ 

N{Y -y,n,ri;b) = 0oe^^ (4.49) 

Comparing Eq. (j4.49p with the functional of Eq. (j4.45p and recalling that ]\f^^^^ exp[^z] 
for ri approaching 1/Qs, we obtain that in the region of long distances, but for ri close to 
1/Qs, that the generating functional for scattering amplitude has the form: 

Ni,iY-y,n,b;[^{r,)]) = [ rfV, 7(r,) e"^ ^^^^'^ (4.50) 



Using Eq. (|4.50|) and Eq. (j4.48p we can see that in the region of large z (Y — y ^ 
ln(l/r^) ^ l)the generating functional is 

Ni,{Y-y,n,b-[^{r,)]) = 1 - )' (4.51) 

with the functional $([7 (?"«)]) defined as 

$([7(r,)]) = ln(^| tiV,e-^"^(i/'-?) 7(n)) (4.52) 

where z is defined by Eq. ()4.47j) . 

The accuracy of Eq. ()4.51|) is not very high and we cannot guarantee the value of the 
constant in front of the exponent in Eq. ()4.51|) . 

4.1.2 Dipole-dipole scattering amplitude in lancu - Mueller approach: 

To calculate the dipole-dipole scattering amplitude we would like to generalize Eq. ()3.39p . We 
replace 7(rj) by z/, while we choose 7(r^) = 7(rj,r9/i^ (see Eq. fl2.29p ). 

Using these functions 7 we can rewrite Eq. (j2.25p in the form of Eq. (j3.39|) . namely, 

Nir,R,Y;b) = 7^ <f — N (Y - y,r,b;[u]) N (y, R,b;[^{r,,r'^ /u]) (4.53) 
Zm JC2 V 
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where we deform contour Ci to contour C2 (see Fig. [7|). In Eq. (j4.53|) we use r and R instead 
of Ti and r2 denoting the sizes of colhding dipoles. 

Substituting Eq. ()4.5H) in this equation we obtain the resulting expression 

N{r,R,Y;b) = (4.54) 



2 vr i Jc2 y 

In Eq. ()4.54j) we denote the sizes of interacting dipoles by r and R instead of r\ and r2- As we 
have seen in toy model example, the contour of integration is situated between all singularities 
of two amplitudes in Eq. ()4.54|) . Due to this fact the contributions such as 

_e-| ( ^^+*[7{n)] and - e-H 

are equal to zero since they have singularity outside of contour C^- Since the integral over the 
large circle (see dotted line in Fig. [Tj) is equal to zero, Eq. ()4.54|) can be rewritten as follows 



N{r,R,Y;b) = (4.55) 



1 / 9\ r+°° 



1 + - ) / dl e 

71 



We are not sure of the value of the coefficient in front of the integral, but we wrote it explicitly, 
to show that the second term is negative. In Eq. ()4.55|) we extract the variable v from both 
<J>'s and use the notation / = In z/. Taking the integral over / explicitly one obtains 

N{r,R,Y;h) = 1 - e-Te-'i^'^-'R) (4.56) 

where 

z{Y,r,R) = AasY - ln(i?V^^) (4-57) 

Eq. ()4.56|) reproduces the result of the lancu and Mueller paper fT] ( with numerical 
coefficients of our simplified model for the BFKL kernel given by Eq. ()4.42|) ) and it does not 
contain any variable related to the fictional dipoles with rapidity y as well as rapidity y itself. 

The sum of 'fan' diagrams leads to the answer [22] 

Nf,n{r,R,Y;b) = 1 - e" i (4.58) 

which is much closer to unity than the correct answer given by Eq. ()4.56|) . In the case of fixed 
QCD coupling our calculations as well as the lancu and Mueller ones, show that the non-linear 
BK equation is not able to describe the physics of the dipole-dipole interactions, even when 
dipoles have quite different sizes (say i? r). 
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4.2 Running a^: 



It is not theoretically clear how to include the running QCD coupling in this approach, 
mostly because we do not know how to write the BFKL kernel for running QCD coupling. 
However, for short distances in the pQCD phase {r"^ Ql < 1) we know that in the BFKL 
equation depends on the size of the scattered dipole [3I211I2S1- This dependence leads to the 
saturation scale QsiY) oc and to the geometric scaling behaviour versus the variable 

r^Qs (at least within the semi-classical accuracy) El! • 

Inside the saturation region {r"^ Ql ^ 1) we can assume that the running is frozen at 
the saturation scale This assumption appears natural from the point of view that physics 
of saturation is determined by one scale: the saturation momentum [3 El EI- Indeed, it was 
shown in Ref. ^31, that this assumption leads to the geometrical scaling j2S] in the saturation 
region, making our hypothesis self-consistent and providing the natural matching with the 
pQCD region. 

In the framework of this assumption we obtain the same solution given by Eq. ()4.4fi|l but 
with a new variable z, namely 

z = In (Qlrl) = 2as{Ql)iY - y) - Hr^/rl) + 2/3(6) . (4.59) 

In Eq. ()4.59p we use the expression for the saturation scale in the case of running 0:5 obtained 
in Refs. El EH]. We would like to recall that in Refs. El CH] instead of factor 
2 in the second line stands (1/2) dx{'y)/d'j\^=^^^, where 7cr is the solution to the equation 
x(7cr)/(l — 7cr) = — c?x(7)/'^7l7=7cr- However, in DLA the BFKL kernel xil) is so simple that 
7cr = 1/2 and (1/2) dxi'j) / d'j\^=^^^ = 2. In section 4.3 we will discuss the general form of the 
BFKL kernel and all these numerical factors will reappear in our calculations. 

Finally, after performing the same calculation as has been discussed above, we obtain the 
same formula of Eq. ()4.56p but with 

Z{Y,r,R) = 2as{QsiY-y)){Y-y) + 2as{Qsiy))y - HR^r') 

= ^ {^/y^ + Vy) - HRVr') ; (4.60) 

The minimum of Z{Y, r, R) occurs aX y = or aX y = Y . The case y = in our notation 
corresponds to the 'fan' diagrams of Fig. El-a type (see also Fig. 0), if we assume that R > r. 
We indeed assumed this, considering ln(rj/r) > and ln(_R/r^) > when we tried to justify 
the generating functional approach for the 'fan' diagrams. As has been discussed in section 
2.2 (see also j3] ), the selection of the 'fan' diagrams was defined for the dipole sizes less than 

1/Qs. 

Therefore, only the solution for y = can match the amplitude at r < 1/Qs- Finally, our 
calculations support the principal idea of the non-linear equation that in the case of running 
QCD coupling it describes the interaction of the dipoles with different sizes. However, this 
conclusion is based on an additional assumption that the running QCD coupling is frozen on 
the saturation scale. 
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As was argued in Ref . |^ , in the approach given by Eq. (j4.42p , it seems natural to assume 
that as depends on the size of produced dipole. Indeed, in this case we expect that the size 
of the produced dipole (r' or |r — -PI in Eq. ()2.23|) ) is smaller than the size of the scattered 
dipole but larger than 1/Qs inside the saturation domain j22[ l25]. Therefore, we can replace 

/ ^^7^^ ^ ^ / . . ^s{r')— + vr / a^dr-r-l) A. , . (4.61) 

Introducing a new function N{C,,y) = /I" as{^') dC,' N{xi',y) where is the dipole 
scattering amplitude and = ln(l/(r^A^)) and = In (Q^/A^), we can rewrite the Balitsky- 
Kovchegov equation in the following way 

which is a direct generalization of Eq.(2.13) in Ref. [22] • Replacing A^(^, y) by function y) 
using A^ = 1 — e~'^^^'y^ one can see that Eq. (j4.62p reduces to 

dy TT 

Assuming that increases moving inside of the saturation region , one can see that A^ ap- 
proaching 

N yl-(l) (4-64) 

Substituting Eq. ()4.64p in Eq. ()4.63|) and integrating with respect to y using the explicit 



expression for i^s = y y) we obtain 



= '-^^^f^{y (^l(ln,-l) + ln^l^- Inel + 0(0} (4.65) 

Function C(^) in Eq. ()4.65p can be found from the condition that </)(^ = ^s) = Const . The 
final answer has the form ^ 

0(e, y) = ^ {^s (HUO - 1) + e) (4.66) 

One can see that this solution does not show the geometric scaling behaviour. However, if 
^ is close to Eq. ()4.66p degenerates into geometrical scaling behaviour within accuracy 
ln(Q^ r^)/ln(Qg/A^). Therefore, this solution is not worse that the one discussed before with 
the exact geometrical scaling behaviour. 

^This solution is very similar to the solution obtained in Ref.j25_ but it has an extra factor as{£,)- It should 
be stressed that only with this factor this solution could be matched with the geometrical scaling solution in 
the pQCD region (r^Q^ < 1 ). 
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Repeating the same calculation as in section 4.1 but with the solution given by Eq. (j4.65|) 
we obtain the behaviour 

N{r,R,Y;b) = 1 - e-'^^^-?^''^)-'^^*''^) (4.67) 
One can see that minimum of the sum of 0's occurs at 

esiVrmn) = ^^(r) (^^^ (4.68) 

where we used obvious notation C,s{Y) = \n{Ql{Y)/A),C,{r) = ln(l/(r^A^)) and so on. 
Assuming that as{r) <^ as{.R) we obtain the following behaviour for the scattering amplitude 
in lancu - Mueller approach: 

^{Y-y,r) + <Piy,R) <P{Y,r) + ^ (-C'(2/™„) +eiR)) , (4.69) 

where ymin is defined by Eq. fl4.68|) . 

One can see that ymin depends on Y and, therefore, the asymptotic behaviour of the 
scattering amplitude cannot be calculated just using the Balitsky-Kovchegov equation on the 
contrary to our example with frozen as- On the other hand, for as{r^ -C O-siR) Vmin is 
almost constant and the influence of the target is rather small. 

These two examples illustrate the importance of understanding the argument in QCD the 
running coupling constant in the framework of the non-linear equation. In our simple model 
with Eq. (j4.42|) for the BFKL kernel the second example looks more reliable. 

4.3 High energy scattering amplitude with the general BFKL ker- 
nel 



The estimates discussed above are based on the solution to the non-linear Balitsky-Kovchegov 
equation given in Ref. |221- In this section we are going to show that this solution cannot give 
a correct high energy asymptotic behaviour. In this section we use the full BFKL kernel and 
solve the Balitsky-Kovchegov equation in the saturation region. 



4.3.1 The solution to the non-linear equation in the saturation region: 

To search for the solution in the saturation region, we use several ideas and technical methods 
that have been discussed in Refs. [221 HSI- 

1. It is useful to consider the non-linear equation in a mixed representation, fixing impact 
parameter 6, but introducing the transverse momenta as conjugated variables to the dipole 
sizes. The relations between these two representations are given by the following equations 

[1311211 

N(r,y;b) = rM kdk Mkr) N(k,y;b) ; (4.70) 
Jo 

N{k,y;b) = / - J,{kr) N{r,y;b) ; (4.71) 
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2. In this representation the non-hnear equation reduces to the form: 

^^g^^'^^ = as {xm))Nik,y;b) - N\k,yM) (4.72) 
where x{l{C) is an operator defined as 

7(0 = 1 + ^ (4-73) 

where ^ = ln(A;^ k'^b^), and k and k' are the conjugated variables to the dipole sizes of the 
projectile and the target, b is an impact parameter which we assumed to be large in Eq. ()4.72|) . 

3. We expect there to be geometrical scaling behaviour of the scattering amplitude in the 
saturation domain. It means that N{k,y] b) is a function of the single variable 

z = HQl{y,b)/A') - m = asj^^iy - Vo) - ; (4-74) 

J- O'er 



where 7cr ~ 0.37 is a solution of the equation 01211 



(4.75) 



djcr 1 - 7cr 

4. Introducing a function 0(z) we are looking for the solution of the equation in the form 

N{z) = ^J'dz'[l - e-<^(^')) ; (4.76) 

5. We assume that function is a smooth function, such that (pzz ^ <Pz (pz where we 
denote 0^ = d(f)/dz and (pzz = d"^ ip / {dzY . This property allows us to rewrite 

e-^^'^ = {-(t)zT e-^^'^ (4.77) 



{dzY 

Substituting in Eq. ()4.72|) N in the form of Eq. ()4.76|) and replacing ?/ by z we obtain 

as = «5 ( [/ X(l - /) - 1] N{z) + N{z) e-^ ) (4.78) 

after taking the derivative with respect to z on both sides of Eq. ()4.72|) . / denotes f = d/dz 
in Eq. ()4.78|) . An important property of function /x(l ~ /) ~ is the fact that at small / 
it has the expansion that starts ^ from Using Eq. ()4.77|1 one can see that the first term on 

■^It should be stressed that the simphfied model for x(l ~ f) — f(^i-^f) does not have this property. This is 
an explanation why in Ref. "22" where this model was used, the solution was missed as we will discuss below. 
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r.h.s. of Eq. (j4.78|) is proportional to e Canceling e ''''on both sides of Eq. (j4.78p and 
once more taking the derivative with respect to z we reduce Eq. ()4.78|) to the form: 



(4.79) 



L(0.) = ^-^^^ -/-^ - S (4.80) 

The answer in the simple model given by Eq. ()4.42j) obtained in Ref. [22] corresponds to the 
simplified version of x(l ~ /) = !//• In this case the second term in Eq. ()4.79j) vanishes and 
the solution for large z is 

For the full BFKL kernel function is negligibly small for small (j)z- However, 

being small for 0^ < 1 has a singularity at — 1 namely 



dU 



d( 



^zf 



for 



The solution given by Eq. 1)4.8111 is controversial since it was obtained assuming that ^^^^"-^ 
is small but it leads to large 0^ where this term is essential. We can conclude that this solution 
is valid only in the simplified model with the kernel given by Eq. ()4.42|) . 

For the full BFKL kernel we can try an opposite approximation assuming that (j){z) is large 
(0 ^ 1) but 0;j is approaching to 1 at z >> 1. In vicinity 0^ = 1 we reduce Eq. ()4.79|) to 

1 d^(h 



(1 - 0,)2 ^dzf 
For large z Eq. ()4.82|1 has a solution 



1. (4.82) 



0(5) = z - \nz (4.83) 

which can be verified by explicit calculations. It should be stressed that 0(z) of Eq. ()4.83|) 
satisfies all conditions of a smooth function that has been used for the derivation of Eq. ()4.79|) . 
We can also check that the l.h.s. of Eq. ()4.79|) is proportional to 1/ z^ and it can be neglected. 

It is more convenient to use the coordinate representation for the scattering amplitude ( see 
Eq. ()4.7()j) and Eq. ()4.71)1 ). We can simplify Eq. ()4.71)1 in the saturation region where N{r, y; b) 
manifests a geometrical scaling behaviour being a function of one variable N{r'^Q'^{y, b)). The 
main contribution in this equation stems from the kinematic region: kr < 1 and r^Q^ > 1. 
Indeed, for kr > 1 the integral of Eq. ()4.7H) is small, due to the oscillating behaviour of 
jQ^kr) while for r^Q^ < 1 the amplitude being in the perturbative QCD domain is rather 
small. Therefore, Eq. ()4.7ip can be rewritten in the form: 

N{k,y;b) = N (z) = ^ J' dz'N{z') (4.84) 
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where z' in the integral is defined by Eq. (j4.47|) . Comparing Eq. (|4.76j) .Eq. (j4.83|) and Eq. (j4.84|) 
one concludes that 

N [r'^Ql^y,}))) = 1 - e-^W + i'^^W (4.85) 

with z{r) is defined by 

z{r) = HQlr') = asj^^ {v - Vo) - Hl/r'A') . (4.86) 

1 - Icr 

We are aware that the majority of experts, including one of us, have thought that solution 
of Eq. ()4.81|) gives a correct asymptotic behaviour for the scattering amplitude. Indeed, at 
first sight, the approach with the kernel of Eq. (j4.42j) is well motivated while the same 
result can be derived without addressing any model for the kernel ( see Ref. |T| and references 
therein). Function has a very transparent physical meaning (see Eq. ()4.46p ): exp""^ is a 
probability for a dipole to pass the target without any inelastic interaction. In the solution 
given by Eq. ()4.81|) . exp""^ is the probability that no extra gluon could be emitted (see Fig. |H1- 
a). However, this solution does not suppress the elastic scattering shown in Fig. |Hl-b. This 
diagrams can easily be calculated 

« ^ /j f (i) - ^ ("^) 

where /Ql describes the emission of the 'soft' gluon with 1/r < q < Qs from the dipole with 
the typical size 1/Qs- Since N^i is the amplitude at fixed impact parameter we divided the 
cross section by the area of the scattering dipole. The idea that approaching to the black disc 
will be described by the Born approximation in the saturation region has been discussed j2Zl 
but we found that this asymptotic behaviour follows from the Balitsky-Kovchegov non-linear 
equation. 

A kind of argument against could be the result of the numerical simulations for dipole- 
dipole scattering by Salam in Ref. j2] which supports Eq. (j4.8H) . in spite of the fact that 
the general BFKL kernel has been used. In Fig. IHl we compare two solutions (Eq. (j4.8ip and 
Eq. ()4.85p ) with the numerical solution of the Balitsky-Kovchegov equation given in Ref. |28j. 
We fit the numerical solutions by two function: 

N{z) = 1 - Cl ; (4.88) 

N{z) = \ - Cl e-(--'"(-+^l)) ; (4.89) 

One can see that both formulae could describe the solution. Therefore, we cannot use the 
numerical solution as an argument in favour of any solution. It is interesting to mention that 
the value of in Eq. (|4.88|) (C2 = 0.11) is closer to the theoretical value in Eq. (j4.8ip than 
it is found by Salam. 



4.3.2 High energy amplitude in the saturation region: 
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Qs > qj >qi^i >i/r 



Qg > q > 1/r 



a) b) 

Figure 8: The gluon emission which is suppressed in the solution given by Eq. (|4.8ip (Fig.|S]-a); 
and the elastic rescattering (Fig. |Hl-b) which is taken into account in the solution by Eq. (j4.85|) 



Fixed a 5 : 

We can calculate the scattering amplitude in the saturation region using the lancu - Mueller 
factorization in the form of Eq. ()4.5H|1 but with the contour Ci for integration over u. For 
N{Y — y,r,b; [z/]) we use Eq. fj4.85|) with z replaced by 2; + 2$[7(rj)] where $ is defined in 
Eq. ()4.52|) and z = In — ln(l/r^). As in the case of the double log kernel we choose 
7(ri) = z/ and 7(r9 = ^{ri,r'^/u. 

We can easily take the integral of Eq. ()4.53|) if we neglect the u dependence in In z. In this 
case the product N(Y — y,r,b[i'] N{y, R,b; ['j{ri,r'j)/i']) does not depend on u. This simple 



Figure 9: The fit of the numer- 
ical solution to the Balitsky- 
Kovchegov equation 28j using 
Eq. and Eq. The 

dashed line corresponds to 
Eq. (jCTjl with Cf = 0.69 and 
C| = —1.05 while the dashed- 
dotted line describes Eq. ()4.88|) 
with Cl = 0.25 and = 0.11 . 
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observation leads to the answer for the scattering amphtude 

N{r,R,Y;b) = 1 - e" ^(^'^'^^^^ in{.{y,r,R)/2) (490) 

where z{Y, r, R) is given by Eq. ()4.57|1 . 

As in the case of the simphfied kernel of Eq. ()4.42j) , the main contribution stems from the 
set of enhanced diagrams. 

Running as '• 

Calculating the diagram of Fig. |Hl-b for the running QCD coupling we obtain the answer 
^(y,r) = Hr'Q-i) - ,„ (,„ (^-^ip^)^)) (4.91) 

with Qs oc e^/^. We can also verify that Eq. (j4.91|) is the solution to the generahzation of the 
equation of Eq. (j4.72j) with the prescription for running as given in Ref. j25j . 

This change leads to remarkable consequences. Indeed, since the minimal value of Z(Y — 
y,r) + Z{y,R), where Z{y,R) ln(Q^(?/) i?^), occurs at?/ = OorF — ?/ = the resulting 
amplitude in the saturation region is 

Nir,R,Y;b) = 1 - e" (4.92) 

For simplicity, we did not write in Eq. (|4.92j) the lnln(. . . in Eq. (j4.92p . 

Therefore, the asymptotic behaviour of the scattering amplitude for running 0:5 is deter- 
mined by the non- linear equation Oil], which is able to describe the colour glass condensate 
phase of QCD [Sj in this case. 

5 Conclusions 

In this paper we showed that the lancu- Mueller factorization is closely related to a sum of a 
particular set of enhanced diagrams (see Fig. ^ in the Reggeon -like diagram technique based 
on the BFKL Pomeron. This set of diagrams leads to the high energy scattering amplitude in 
wide range of energy and/or Bjorken x, namely 

1 

Y = ln(l/x) < — In . (5.93) 

as as 

We found a simple formula of Eq. ()4.53|) which sums these diagrams and manifests the 
lancu-Mueller factorization. Using this formula we confirm the lancu-Mueller result, that 
at fixed QCD coupling, the high energy asymptotic behaviour is determined by the enhanced 
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diagrams, but not by the 'fan' diagrams that are summed by the means of Bahtsky-Kovchegov 
non-hnear equation. We showed that the dipole -dipole amphtude in the saturation region 
can be written as 



However, for the running QCD couphng the main contribution in the saturation region stems 
from the 'fan' diagrams of Fig. |21-a (or Fig. E)) -type. This new result allows us to penetrate 
the saturation region, summing the 'fan' diagrams of more general type (see Fig. |21-b for the 
first such diagram) . The first suggestion of how to sum such diagrams was proposed in Ref . 
jl4j . and the present paper confirms our hope that the general 'fan' diagrams will lead to the 
estimates of high energy asymptotic behaviour for running QCD coupling in the saturation 
region. 

Our final result for the asymptotic behaviour of the scattering amplitude is quite different 
from that in the lancu - Mueller paper even for fixed QCD coupling case. This difference 
stems from the new solution to the non-linear equation found in this paper, which includes 
the general BFKL kernel as opposed to the simplified version of the kernel ( see Eq. ()4.42j) ) 
which has been used in Ref. ||2'2|j . 

Our approach to the lancu-Mueller factorization is even closer to the Monte Carlo simula- 
tion than the original papers jT], and the semi-analytical estimates suggested in Ref. We 
hope that the simple formula suggested in this paper will give an impetus for searching for a 
new approach to find the high energy amplitude in the saturation region. We believe that we 
have demonstrated in this paper, that the lancu-Mueller factorization provides a method to 
take the essential fluctuations in the partonic wave function of colliding dipoles into account 
and shows the practical method of calculation beyond of the non-linear equation. It also allows 
us to find the region of applicability of the non-linear equation. 
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